function s = numbersteps(E,N,share,alpha,gamma,psi)

A = getA(E,N,alpha,gamma);
s = log((A.*share).^(1-alpha))./log(psi);

end
